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Mathematical  Modeling  of  the  Temperature 
Rise  in  a  Thin  Cell  Culture  Exposed  to  High 
Frequency  Electromagnetic  Irradiation 


1  Introduction 

This  report  describes  a  mathematical  model  to  predict  the  temperature  increase 
in  a  thin  layer  of  tissue  culture  exposed  to  RF  energy.  This  effort  is  motivated  by 
an  experiment  conducted  by  Dr.  Nancy  Millenbaugh  in  the  Air  Force  Research 
Laboratory  at  Brooks  City  Base,  Texas.  To  put  this  mathematical  model  in  the 
proper  perspective,  we  will  describe  briefly  the  biology  involved.  Heat  stroke 
(sudden  fall  in  arterial  blood  pressure)  can  occur  under  severe  hyperthermia, 
leading  possibly  to  circulatory  shock  (tissue  hypo  perfusion).  The  exact  mech¬ 
anism  is  still  in  question.  Research  suggests  that  the  loss  of  vasoconstrictor 
tone  and  peripheral  pooling  of  blood  may  be  involved  in  heat-induced  circu¬ 
latory  shock.  Similar  phenomena  were  observed  when  heat  stress  is  replaced 
by  millimeter  wave  (MMW)  heating.  (The  MMW  band  is  from  30  to  300  Hz, 
corresponding  to  a  wavelength  of  10  mm  to  1  mm.)  However,  Freni,  et  al.  [4] 
suggest  that  environmental  heating  and  MMW  heating  are  not  the  same.  For 
one  thing,  core  temperatures  under  MMW  stress  do  not  rise  as  much  as  that 
under  environmental  heat  stress,  when  circulatory  failure  begins  to  occur.  For 
another,  MMW  exposure  can  produce  larger  and  more  rapid  increase  in  subcu¬ 
taneous  temperature  than  environmental  heating.  They  suggest  that  cutaneous 
thermoreceptor  may  be  responsible  for  the  initiation  of  circulatory  failure  dur¬ 
ing  such  extreme  peripheral  heating.  For  the  human  body,  most  of  the  MMW 
absorption  is  in  the  region  of  the  cutaneous  thermal  receptors  (0.1  mm  -  1.0 
mm)  [5]. 

To  discover  possible  mechanisms  involved  in  MMW-heating  in  biological  en¬ 
tities,  an  experiment  at  Brooks  has  been  conducted  to  expose  cells  growing  in 
culture  flask  by  35  GHz  MMW.  By  submitting  the  exposed  cells  for  differen¬ 
tial  gene  expression  analysis,  the  researchers  hope  to  find  clues  to  mechanisms 
underlying  MMW-heating. 

2  Some  Experimental  Results 

The  flask  containing  the  cell  culture  is  shown  in  Figure  1.  A  transmitter  gen¬ 
erating  a  35  GHz  plane  wave  is  coming  from  the  bottom  of  the  flask  and  is 
polarized  along  its  long  axis.  Using  an  infrared  camera,  the  temperature  in  the 
cell  culture  can  be  monitored  and  captured.  An  example  of  the  temperature 
excursion  in  a  cell  culture  is  shown  in  Figure  2.  Another  example  is  shown  in 
Figure  3.  However,  here  a  jelatin-like  substance  (whose  heating  characteristics 
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are  expected  to  be  similar  to  that  of  the  actual  targeted  cell  culture)  is  used 
instead.  Figure  4  shows  the  temperature  rise  in  the  same  but  at  two  selected 
cross-sections  of  the  flask.  Researchers  would  like  to  sample  cells  for  later  ge¬ 
netic  analysis  at  locations  where  the  temperature  rise  is  insignificant  in  order 
to  minimize  temperature  effects.  This  would  afford  them  a  better  chance  of 
discovering  possible  non-thermal  effects  of  EM  radiation.  The  question  is  there¬ 
fore  where  is  the  best  place  to  sample  the  cells?  The  question  is  not  easy  to 
answer  because  the  observed  temperature  excursions  seem  to  be  sensitive  to  the 
experimental  setup  as  the  observed  temperature  excursions  are  often  different 
from  “run”  to  “run”.  This  is  especially  true  in  the  earlier  experiments  (shown, 
for  example,  in  Figure  2).  This  could  be  due  a  small  change  in  the  flask  position 
in  the  radiation  field  or  in  other  unknown  factors. 


Figure  1:  Cut-away  flask  with  cell  culture 

The  goal  of  the  mathematical  modeling  effort  described  here  is  to  try  to 
understand  what  has  been  observed  experimentally.  In  particular,  what  tem¬ 
perature  excursions  are  to  be  expected  in  the  cell  culture  or  its  jelatin-like  simu¬ 
lant,  if  the  flask  is  indeed  incident  by  a  (35  GHz)  plane  wave  orthogonally.  Our 
approach  is  to  model  the  experimental  setup  mathematically.  We  first  calculate 
the  specific  absorption  rate  (SAR)  and  then  use  it  in  turn  to  calculate  the  tem¬ 
perature  excursion  in  the  cell  culture  or  its  jelatin-like  simulant.  (For  simplicity, 
we  will  not  distinguish  between  cell  culture  and  its  jelatin-like  simulant  in  the 
following.) 
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Figure  2:  Temperature  Excursion  at  approximately  4-second  intervals.  Black 
is  “cold”  and  white  is  “hot”. 
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Figure  3:  Temperature  Excursion  at  approximately  3.5-seconds  intervals  in 
jelatin-like  substance.  Black  is  “cold”  and  white  is  “hot”. 


Figure  4:  Temperature  rise  in  a  jelatin-like  substance  in  selected  cross-sections. 


3  SAR  Calculation 

Both  the  thinness  of  the  cell  layer  and  the  high  frequency  of  the  incident  wave 
require  us  to  use  small  spatial  resolution  (Ax  =  150E-06  m).  This  in  turn 
requires  us  to  use  very  small  time  steps  (approximately  290  fs).  To  simplify  the 
calculation  here  and  that  of  the  temperature  later  on,  we  model  the  cell  culture 
as  being  housed  in  a  shallow  rectangular  flask  (Fig  5). 


Figure  5:  Cell  culture  in  a  shallow  flask. 
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We  used  a  commercial  finite  difference  time  domain  code  (XFDTD)  to  cal¬ 
culate  the  SAR  in  the  cell  culture  when  it  is  incident  by  a  35  GHz  plane  wave 
with  an  amplitude  of  86  V/m.  The  above  mentioned  spatial  resolution  results 
in  approximately  14  million  Yee  cells  (Nx  x  Ny  x  Nz  =  651  x  593  x  36).  Fortu¬ 
nately,  steady  state  was  achieved  in  a  relatively  short  time  (approximately  0.4 
ns,  or  1,400  time  steps). 

The  electric  parameters  used  for  the  flask  and  the  cell  culture  are  shown  in 
Table  1.  The  average  SAR  (W/m3),  denoted  henceforth  by  Sv,  at  a  point  r  in 
the  flask  or  in  the  cell  culture  is  calculated  using  the  formula 

5v(r)=cr(r)||E(r)||2ms  (1) 

where  ||E(r) ||rms  is  the  root  mean  square  value  of  the  electric  field  (V/m)  and 
<j(r)  the  conductivity  (S/m)  at  the  point  r.  Typical  Sv  values  are  shown  in 
Figures  6-7.  Figure  6  shows  the  values  of  Sv  in  a  layer  that  includes  both  the 
flask  and  the  cell  culture.  Figure  7  only  shows  those  in  the  cell  culture.  It  is 
seen  that  the  highest  Sv  values  tend  to  be  in  the  interface  between  the  two 
media. 

In  most  EM  dosimetry  studies  the  estimation  of  Sv  is  the  endpoint.  However, 
for  our  investigation  here,  this  is  only  the  beginning,  as  we  are  interested  in 
knowing  the  temperature  distribution  throughout  the  cell  culture  beyond  the 
initial  exposure.  (It  is  well-known  that  the  initial  temperature,  for  example, 
during  the  first  few  nanoseconds,  is  directly  related  to  the  average  SAR.) 


x  10 


y  0  0 


Figure  6:  Typical  SAR  in  a  layer  including  the  flask  and  the  cell  culture. 
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Figure  7:  Typical  SAR  in  a  layer  inside  the  cell  culture  only. 


4  Temperature  Rise  Calculation 

4.1  3D  Diffusion  Equation 

To  calculate  the  temperature  rise,  we  start  with  the  general  3D  heat  conduction 
equation 


-V  -q  +  Q(x,y,z,  t) 
q(x,y,z,t) 


-kVT(x,y,z,t) 


(2) 

(3) 


Or,  simply, 


1  9T  1  ^ 

— 3-  =  V  T  +  -Q 
a  dt  k 


Here  q  =  (q1,q2,q3)  is  the  heat  flux,  wherein  each  component  g,  (i  =  1,2,3) 
is  measured  in  W/m2,  k  is  the  thermal  conductivity  (W/m/K),  p  the  density 
(kg/m2),  and  cp  the  specific  heat  (J/K/kg).  The  parameter  a  is  the  thermal 
diffusivity  (m2/s)  and  is  defined  by 


_  k_ 
a~  Pcp 

Finally,  Q  is  the  heat  source  (W/m3)  for  the  problem. 


4.2  The  3D  Code 


The  3D  code  we  used  is  a  modified  version  of  a  more  general  diffusion  code 
obtained  from  J.  Zhang  and  J.  J.  Zhao  at  the  University  of  Kentucky  [25]. 
They  use  it  to  solve  a  special  3D  microscale  heat  transport  problem  described 
by  the  equations  (assuming  the  scale  in  the  z-direction  is  small  in  the  order  of 

0.1  H) 


—V  •  q  +  Q 

dT 

~  pCp  dt  ’ 

(4) 

?!  (*,»,«!*) 

=  ~kTw(x,y,z,t), 

(5) 

q2(x,y,z,t) 

=  -kTy(x,y,z,t), 

(6) 

Q3(x,y,z,t  +  Tq) 

=  ~kTz(x,y,  z,t  +  tt) 

(7) 

where  Q  a  heat  source,  which  in  our  case  is  a  |  E  |2  . 

The  only  differences  between  Equation  (3),  the  Fourier’s  Law  of  heat  conduc¬ 
tion,  and  Equation  (7)  are  the  time  delays  rq  and  rT.  The  reason  for  these  time 
delays  (also  called  relaxation  times)  is  basically  to  account  for  the  finite  speed 
of  heat  propagation  [18].  As  the  relaxation  times  are  related  to  the  reciprocals 
of  speeds  of  heat  propagation,  the  classical  model,  corresponding  to  the  case  in 
which  the  "relaxation”  times  rq  and  rT  vanish,  implies  an  infinite  speed  of  heat 
propagation.  The  insistent  of  finite  propagation  speed  is  especially  important 
in  microscale  heating,  though  it  may  be  important  in  our  case  also. 

Equations  (5-7)  are  clearly  related  to  the  dual-phase-lag  model  of  Tzou  [18]: 

q(x,y,z,t  +  rg)  = -kVT(x,y,z,t  +  rT)  (8) 

which  in  turn  is  related  to  many  existing  microscale  heat  transfer.  When  rT 
vanishes  in  the  dual-phase-lag  model,  Equation  (8),  one  basically  gets  the  well- 
known  Cattaneo  and  Vernotte  (CV)  equation  [2,  20] 


q  +  r^=-kVT  (9) 

which  together  with  Equation  (4)  leads  to  the  hyperbolic  heat  conduction  equa¬ 
tion 

dT  d2T 

- h  t - =  aV2T 

dt^  dt2 


This  equation  is  also  known  as  the  thermal  wave  model  or  the  telegraph  equation 
[9].  There  are  some  experimental  observations  [14,  10]  that  resemble  hyperbolic 
heat  conduction  in  living  cells  as  well,  although  the  basic  mechanism  is  still 
controversial  [22]. 

In  this  study,  we  have  restricted  ourselves  to  the  simple  case  in  which  rq  — 
0  and  rT  =  0,  and  the  possibility  and  the  significance  of  hyperbolic  heating  in 
high  frequency  irradiation  will  be  investigated  in  the  future. 

The  code  uses  the  Crank-Nicholson  method  to  discretize  the  differential 
equation.  To  solve  the  resulting  large  discrete  system  efficiently,  the  code  uses  a 
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preconditioned  conjugate  gradient  method.  The  preconditioner  uses  the  incom¬ 
plete  Cholesky  factorization  of  the  system  matrix  which  is  positive  definite  and 
symmetric.  All  of  these  computations  are  done  using  a  sparse  representation  of 
the  system  matrix. 

4.3  Boundary  Conditions 

The  solution  depends  on  the  external  boundary  conditions  as  well  as  internal 
boundary  conditions.  When  the  Dirichlet  boundary  condition  is  applied  by  set¬ 
ting  the  boundary  temperature  to  the  constant  external  temperature,  a  steady 
state  is  achieved,  as  an  appropriate  amount  of  heat  is  continuously  being  leaked 
out  to  balance  out  the  absorbed  power.  This  is  contrary  to  what  is  observed  in 
the  experiment.  As  the  original  3D  code  handles  only  Dirichlet  boundary  con¬ 
ditions,  we  are  led  to  implement  a  more  realistic  mixed  (“radiation”)  boundary 
condition  [1,  Carslaw  and  Jaeger]  of  the  form 

‘I  +  MT-r,)-° 

where  k  is  again  the  (internal)  thermal  conductivity,  h  the  “external  conductiv¬ 
ity”  (W/m2/K),  and  Tt  the  external  temperature. 

4.4  Internal  Boundary  Conditions 

Our  problem  is  non-homogenous  consisting  of  the  flask  and  the  cell  culture.  As 
the  original  code  only  treats  the  homogeneous  case,  we  extended  the  code  to 
handle  simple  internal  boundary  conditions  (on  the  surface  separating  the  flask 
and  the  culture  cell)  of  the  form 


T  =  T 
x\  x2 


and 


where  ^  is  the  differentiation  along  the  normal  to  the  surface  separating  the 
flask  and  the  culture  cell.  These  new  conditions  destroy  the  symmetry  of  the 
system  matrix  and  required  us  to  implement  a  less  efficient  sparse  matrix  solver. 


4.5  Source 

It  is  natural  as  well  as  tempting  to  use  the  instantaneous  SAR  as  the  source: 


Q{x,y,z,t)  =  <r|iE(x,2/,z,<)||2 

However,  it  soon  becomes  clear  that  this  is  highly  impractical  due  to  the  high 
frequency  (35  GHz)  involved.  If  we  were  to  make  full  use  of  the  instantaneous 
SAR,  we  would  need  to  use  a  time  step  on  the  order  of  a  few  picoseconds. 
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This  becomes  impractical  to  solve  (especially  on  a  PC),  since  each  step  involves 
solving  a  system  of  approximately  a  quarter  of  a  million  unknowns,  even  when 
we  use  a  moderately  refined  spatial  discretization.  The  second  best  choice  is 
therefore  to  use  the  time-averaged  SAR  defined  in  Equation  (1)  as  the  source, 
i.e., 

Q(x,y,z,t)  =  Sv(x,y ,  z) 


4.6  Parameters 

Since  it  is  difficult  to  find  the  exact  heat-related  parameter  values  for  our  prob¬ 
lem,  we  have  “averaged”  parameter  values  extracted  from  the  following  sources: 

•  Lu,  et  al.  [12]  have  parameters  for  different  human  tissues  (eye,  bone, 
brain,  skin  and  CSF)  at  1  GHz.  The  values  (averaged)  are  (er  =  49.76,  o  = 
1.24  S/m,  and  k  =  0.48  W/m-K) 

•  Van  Leeuwen,  et  al  [19]  have  the  following  parameters  (averaged)  for  brain 
tissues:  er  =  54.89,  a  =  1.22  S/m, k  =  0.48  W/m-K, cp  =  3529  J/kg-K,  and  p 
1014  kg/m. 

•  Wainwright  [21]  has  parameters  for  various  human  tissues.  The  averaged 
values  are:  p  =  1125  kg/m3,  k  =  0.42  W/m-K  ,  cp  =  3210  J/kg-K, 

•  Hirata,  et  al.  [6]  have  values  for  human  eye  tissues.  The  averaged  values 
are:  cp  =  3838  J/kg-K,  and  k  =  0.54  W/m-K 

•  Hurt  [7]  recommends  the  following  electric  parameters  at  35  GHz:  a  = 
60S/m,  and  er  =  22  for  cell  culture  and  a  =  0  S/m,  and  err  =  2.4  for  the 
flask. 

•  For  the  possible  flask’s  electric  conductivity,  Stratton  [17]  has  a  =  1.0JS?  — 

08  S/m  for  celluloid. 

•  From  [8]  we  have  the  following  heat  parameters  for  flask  (semiconductor), 
cp  =  250  J/kg-K  and  p  =  5000  kg/m3 

•  From  [16]we  have  heat  parameters  for  glass:  k  —  0.836  W/m-K 

•  From  Internet  Web  Site  MatWeb.com,  we  have,  for  Polystyrene,  p  = 
1050  kg/m3,  cp  =  1800  J/kg-K,  and  k  —  0.14  W/m-K 

The  initial  values  of  the  model  parameters  are  summarized  Table  1. 


9 


Table  1.  Model  Parameters 


Parameters 

Cell 

Flask 

er 

22 

2.4 

a  (S/m) 

60 

5.74E-04 

k  (W/m-K) 

0.48 

0.14 

P  (kg/m3) 

1070 

1050 

c,  (J/kg-K) 

3526 

1800 

Derived  Parameters 

1.27E-07 

7.4E-08 

2.085 

7.14 

4.7  Initial  Conditions 

Using  the  parameters  in  Table  1  and  an  assumed  constant  initial  (room)  tem¬ 
perature  throughout,  we  were,  unfortunately,  unable  to  produce  temperature 
profiles  that  match  (qualitatively)  the  experimental  observations.  The  same  is 
true  when  we  experiment  with  (many)  different  sets  of  parameter  values.  This 
led  us  to  think  that  the  assumption  of  constant  initial  temperature  may  not  be 
valid.  Consequently,  instead  of  using  a  constant  initial  temperature,  we  use  an 
initial  temperature  profile  that  is  an  approximation  of  the  observed  (see  Figure 
8).  The  difference  between  the  actual  and  its  approximation  is  shown  in  Figure 
9. 


4.8  Final  Parameter  Adjustments 

As  we  cannot  be  certain  if  the  parameter  values  for  the  jelatin-like  substance 
used  in  the  experiment  of  interest  (cf.  Figures  3  and  4)  are  similar  to  that  for  the 
estimated  cell  culture  values  shown  in  Table  1,  we  allowed  ourselves  the  freedom 
to  adjust  these  parameter  values  to  match  the  experimental  results  qualitatively. 
The  most  significant  change  is  the  thermal  conductivity.  The  numerical  results 
we  will  show  in  the  following  used  a  thermal  conductivity  of  24  (W/m/K)  for 
the  jelatin-like  substance  instead  of  0.48  (W/m/K)  estimated  for  cell  cultures 
in  Table  1. 
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itial  temperature  rise  (at  t— 0  sec)  and  the  quadratic 


Absolute  difference:  fit  -  observed 


Figure  9:  Fit  error:  surface  and  projected  views.  Dark  is  low  and  bright  is 

high. 
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5  Simulation  Results 


Simulation  results  are  shown  in  Figures  10-13.  Figure  10  shows  the  temperature 
profiles  at  different  times  along  a  middle  x-section  (bottom  figure)  and  a  middle 
y-section  (top  figure).  Qualitatively,  this  compares  reasonably  well  with  the 
experiment  results  shown  in  Figure  4.  Recall  the  experiment  uses  a  flask  (cf. 
Figure  1),  whereas  the  simulation  uses  a  rectangular  box  (cf.  Figure  5). 

Figure  11  shows  the  temperature  profile  in  a  layer  near  the  top  at  four 
different  times. 

Figures  12  and  13  show  two  different  views  of  temperature  profiles  in  the 
entire  jelatin-like  substance  after  10  seconds.  Figure  12  is  drawn  to  scale. 


020620-08:  Temperature  rise  along  x;  delt^  =  1  s 


Figure  10:  Simulated  temperature  rise  in  selected  cross-sections. 
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Figure  11:  Simulated  temperature  distribution  in  a  layer  near  the  top  (z=24) 

at  4  different  times. 


Figure  12:  Simulated  temperature  distribution  in  the  whole  cell  culture  after 


Figure  13:  Simulated  temperature  distribution  in  the  whole  cell  culture  after 

11  sec.  (Not  to  scale.) 

6  Conclusion/Discussion 

We  used  a  combination  of  a  3D  electromagnetic  code  and  a  3D  heat  diffusion 
model  to  study  the  temperature  rise  in  a  thin  layer  of  jelat in-like  substance  that 
is  placed  in  a  flask  while  incident  by  a  35-GHz  electromagnetic  plane  wave.  This 
is  motivated  by  the  need  to  estimate  temperature  rise  in  cell  cultures  in  similar 
experiments. 

We  used  a  finite  difference  time  domain  (FDTD)  code  to  solve  the  elec¬ 
tromagnetic  portion  of  the  problem.  Due  to  the  high  frequency  involved  and 
the  thinness  of  the  cell  culture,  it  took  about  30  hours  to  obtain  the  steady 
state  SAR  on  our  Sun  workstation.  A  faster  code  like  the  VMAX  (Monopole 
Research)  has  solved  a  similar  problem  in  about  6  hours  on  a  similar  platform. 

The  3D  heat  diffusion  model  we  used  is  a  special  case  of  the  Pennes  bioheat 
equation  [15] 

BT  2 

pcp^  =  kV2T  +  cb(T  -Tb)  +  cr  \  E  \ 

often  used  to  study  temperature  rise  in  human  exposed  to  an  electromagnetic 
field.  The  Pennes  bioheat  model  and  other  more  elaborate  models  [23]  take  into 


account,  among  other  things,  the  role  of  blood  flow  in  the  tissues.  The  special 
case  we  have  here  does  not  have  to  deal  with  this  aspect. 

The  3D  diffusion  equation  solver  is  normally  quite  efficient  for  the  simple 
homogeneous  problem,  but  less  so  for  the  inhomogeneous  case.  In  the  latter, 
one  has  to  deal  with  an  internal  boundary  (here  between  the  flask  and  the  cell 
culture).  Dai  and  Nassar  [3]  used  a  “domain  decomposition”  method  to  solve  the 
problem  by  basically  solving  two  separate  problems,  one  on  each  domain,  and 
used  an  iterative  method  to  match  up  the  solution  on  the  common  boundary. 
The  method  we  employed  is  simpler  but  did  destroy  the  symmetry  and  hence 
the  efficiency.  When  the  domain  of  computation  is  irregular,  then  a  boundary 
element  type  method  may  be  better  [13]. 

There  is  also  a  basic  question  on  what  the  infrared  camera  is  capturing.  The 
color  at  each  point  in  each  picture  frame  obviously  correlates  to  some  temper¬ 
ature.  The  question  is  how  does  a  3D  temperature  distribution  transform  to 
the  2D  distribution  recorded  by  the  camera?  What  type  of  integration  does  one 
need  to  perform  on  the  3D  temperature  distribution  to  obtain  the  2D  distribu¬ 
tion  captured  by  the  camera?  Or,  does  the  infrared  camera  only  capture  only 
the  temperature  distribution  on  the  surface? 

In  a  series  of  papers  J.  Liu  et  al.  [11]  proposed  a  skin  burn  model  that 
is  different  from  the  traditional  Pennes’  bioheat  equation.  It  includes  thermal 
waves  due  to  ’’delay”  or  inhomogeneity  in  skin.  The  model,  called  the  thermal 
wave  model  for  bioheat  transfer(TWBMT),  yields  results  quite  different  from 
that  using  the  Pennes’s  formulation,  especially  for  a  ’’flash  fire”  situation  (high 
flux  heating  with  extremely  short  duration).  This  suggests  the  traditional  heat 
equation  may  not  be  an  appropriate  model  for  MMW  stress. 

The  energy  equation  (2)  is  valid  for  rigid  bodies  only  [9].  It  is  implicitly  in 
the  Pennes’s  bioheat  equation  also.  As  most  biological  organs  cannot  strictly 
be  considered  as  rigid  bodies,  studies  using  the  Pennes’  bioheat  equation  as  well 
as  the  present  study  should  be  re-examined  using  a  more  appropriate  energy 
equation. 

In  the  attempt  to  match  the  experimental  results  for  jelatin-like  substance, 
we  had  to  increase  the  values  of  some  parameters  to  over  and  above  those 
expected  for  actual  cell  cultures.  This  increase  in  parameter  values  allows  us 
to  match  the  temperature  increase  observed  in  the  experiment.  Without  this 
increase  in  parameter  values,  i.e.  if  we  had  used  parameter  values  documented 
in  the  literature  for  tissues,  the  temperature  increase  would  have  been  less. 

Finally,  as  far  as  matching  experimental  results  is  concerned,  the  initial  tem¬ 
perature  value  in  the  jelatin-like  substance  is  one  of  the  most  crucial  components 
of  the  calculation.  It  is  still  unclear  why  the  temperature  is  not  uniform  at  the 
beginning  of  the  experiment. 
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